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We propose a one-dimensional quantum Heienberg spin-2 chain, which exhibits two topologically distinct 
valence bond solid states in two different solvable limits. We then construct the phase diagram and study the 
quantum phase transition between these two states using the infinite time evolving block decimation algorithms. 
From the scaling relation between the entanglement entropy and correlation length, we determine that the central 
charge for the underlying critical conformal field theory is c = 2. 
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I. INTRODUCTION 

Recently there have been considerable interests in the inves- 
tigations of topological ordered states, and the quantum phase 
transitions between them. Since topological ordered phases 
usually do not exhibit conventional symmetry breaking, these 
phase transitions naturally can not be described by the conven- 
tional Landau-Ginzburg paradigm. Despite of the lack of local 
order parameters, tremendous progresses have been made in 
characterizing topological ordered states. Properties such as 
ground state degeneracy, quasiparticle statistics, existence of 
edge states, topological entanglement entropy,[l, 2] and en- 
tanglement spectrum[.i] have been proposed and used to dis- 
tinguish different topological ordered states. In contrast, the 
study of topological phase transitions is still in its infancy, and 
progresses are in demand. 

One dimensional quantum spin chains have been a subject 
of interests for many years. It started with the famous Hal- 
dane conjecture, [ ] followed by the Affleck-Kennedy-Li eb- 
Tasaki (AKLT) construction of the valence bond solid (VBS) 
states and their associated parent Hamiltonians.[ i] More- 
over, the SU(2) AKLT model has been generated by intro- 
ducing q-deformed SU(2) group,[6] supersymmetry,[7] and 
higher symmetric groups, such as SU(n),[i>] SP(n),[y], and 
SO(n).[10] More recently a systematic method for construct- 
ing translational invariant VBS state for the general Lie group 
has been proposed. [ I I ] 

For the 5=1 VBS state, there is an appealing physical 
picture where each spin-1 is decomposed into two "virtual" 
spin-l/2's. Across each valence bond, two neighboring vir- 
tual spins pair into a singlet, den Nijs and Rommelse pro- 
posed a nonlocal string order parameter (SOP) which revealed 
a hidden "diluted antiferromagnetic order". [12] Kennedy and 
Tasaki found an unitary transformation that turns the nonlocal 
SOP to a local ferromagnetic order parameter associated with 
a hidden Z2 ® Z2 symmetry. [13] However, it is extremely dif- 
ficult to generalize such a description to the cases of higher 
quantum integer spin chains. 

In this paper, we present a model for the S ~ 2 chain 
which exhibits two distinct VBS states in different param- 



eter regimes. For one of the states, each spin-2 is decom- 
posed into two virtual spin-l's, and for the other it is decom- 
posed into two spin-3/2's. These virtual spins then pair up 
across every nearest neighbor bonds. In the following, we 
shall refer to these two states as VBSi and VBS3/2 states. 
For an open chain, the ground states have 9- and 16-fold de- 
generacies in these two cases, respectively. Hence these two 
VBS states have different topological order Interestingly, 
there is a continuous quantum phase transition between these 
two VBS states. By analyzing the relation between the von 
Neumann entanglement entropy and the spin-spin correlation 
length, [14, IS] we deduce the central charge associated with 
the critical conformal field theory to be two. We further con- 
jecture that the underlying critical field theory may be de- 
scribed by the level-four SU (2) Wess-Zumino-Witten model. 

The paper is organized as follows. In Sec. II, we re- 
view the properties of VBSi and VBS3/2 states. In Sec. 
Ill, the quantum phase transition for the above S ~ 2 spin 
model is explored using the infinite time evolving decimation 
method.[!( ] The entanglement spectrum around the phase 
transitions is studied and the central charges for the underly- 
ing conformal field theory at the critical line are determined. 
A summary is given in Sec. IV. 



n. TWO DISTINCT VBS STATES OF A SPIN-2 CHAIN 

The model Hamiltonian of the spin-2 chain is proposed as 

H = ^[J2P2{i,i + l) + J-iP?.{i,i + l) + PA{i,i + l)] 
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where Prii^i + 1) is the SU(2) symmetric operator that 
projects the spin states associated with sites i and z + 1 into 
the total spin-T multiplet. The coupling constants J2 and J3 
are all positive. In order to make the paper self-contained, we 
first review the VBSi and VBS3/2 states, respectively. 



A. The VBSi state - AKLT state 

To construct this state, we view each spin-2 as a symmetric 
product of two virtual spin-l's. In the VBSi state, two neigh- 
boring virtual spin-1 form a singlet. The direct product of 
the spin-2 multiplets on neighboring sites can be decomposed 
into a direct sum of the total spin 5 = 0,1, 2, 3, 4 multiplets. 
Due to the singlet pairing of the neighboring virtual spin-l's, 
the total spin 5 = 3, 4 multiplets can not be generated. For 
J 2 = 0, the model Hamiltonian only penalizes the 5 = 3, 4 
two-spin states. Hence VBSi state is the unique ground state 
of Eq.(l). In an open chain, the unpaired virtual spin-l's at 
the two ends are free, and they give rise to the 3x3 = 9 fold 
ground state degeneracy. [5, 13] 

In order to write down the ground state wave function, we 
can use the Schwinger boson representation, and the spin-2 
operators can be expressed as 

St = a|fe„ S; = b\a,, SI = {a\a, - b\h)/2, (2) 

with a local constraint a\ai + h\b = 4. Then the S = 2 AKLT 
VBS ground states can be expressed in a simple form, [17] 
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where each 
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creates a singlet bond com- 



posed of two spin- 1/2 between i and i + 1 sites. Furthermore, 
by re-arranging the creation operators in Eq.(3) and combin- 
ing operators with the same site operators together, |\E'aklt) 
can be written in a matrix product state form straightfor- 
wardly. 
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where {^["J} with to = ±2, ±1,0 are 3x3 matrixes. 
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B. The VBS3/2 State - SO(5) symmetric state 

Instead of splitting a spin-2 into two virtual spin-l's, one 
can also split it into two spin-3/2's. In the following, we 
shall first view the VBS3/2 state as the AKLT state of a larger 
symmetry group, which is equivalent to the 50(5) symmetric 
matrix product state in a two-leg electronic ladder[ lS]. Af- 
terwards we will rephrase everything in terms of the physi- 
cal spin SU{2). As pointed out in Ref.[19], one can view 
the ±3/2, ±1/2 states of a spin-3/2 as the four states of the 
spinor representation of SO{5). Similarly one can regard the 
±2, ±1, states of spin-2 as the five-dimensional vector irre- 
ducible representation (IR) of SO (5). Analogous to decom- 
posing a spin-1 vector IR of SU (2) into two virtual spin-l/2's 
spinor IR of SU{2), we can view the vector IR as the symmet- 
ric component of the tensor product of two virtual spinor IR's, 
i.e.. 
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The numerals are the dimensions of the SO{5) IR's. The ten- 
sor product of two 5's on adjacent sites decomposes into 

5®5=1®10®14. (7) 

Comparing expressions of Eq.(6) and Eq.(7), we view Eq. (6) 
as the tensor product of two neighboring virtual spins after 
their respective partners have form SO{5) singlet with other 
virtual spins. Then one can find that SO{5) singlet 1 and the 
antisymmetric 10 appear in the decomposition but the sym- 
metric 14 is absent. Therefore, if H = Pi4{i, i + 1), the 
VBS3/2 state where neighboring virtual 4's pair into SO{5) 
singlet will be the ground state. In an open chain, due to the 
unpaired free 50(5) spinors at two ends, the ground states are 
4 X 4 = 16 fold degenerate. A clear and detailed argument 
of this degeneracy is given is section III(B). Furthermore, the 
projection operator Pi4{i, i + 1) can be expressed in terms of 
the 50(5) generators[ 1 0] 

ni(M-) = ^ E Lfl^f + ^i E LfLff + l. 

l<a<b<5 l<a<b<b 

(8) 

Because the physical spin is SU{2), which is a subgroup of 
50(5), each IR of 50(5) must decompose into an integral 
number of SU (2) multiplets. Thus the 14 discussed above 
must be expressible as the direct sum of SU{2) IR obtained 
by decomposing the direct product of two 5 = 2 multiplets. 
Since the 14-dimensional IR is symmetric upon the exchange 
of site indices, it must only contain even-spin SU{2) multi- 
plets. A simple calculation shows that 14— >5 = 205 = 4. 
Consequently, Pi4(z, i + 1) reduces to Eq.(l) with J3 = 0, 
which is first given by Tu, Zhang, and Xiang.[ 1 0] 

Moreover, the 50(5) generators can be represented by 
L"' = Ea,/? V'lr°^^^/>0, where „ creates a spin-3/2 
fermion with spin index a = ±3/2, ±1/2, F", (a = 1, • • • , 5) 
are the 4-dimensional Dirac F matrices, and F"*" = ^[r°-, F^]. 
By using the above representations, the VBS3/2 state can be 
written as [10] 
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where Ps=2{j) is the spin-quintet projector and 
^afs"^] a'^afsi^j^i fj IS an S0{5) invariant valence bond 
singlet creation operator. TZ is the 5*0(5) invariant matrix. 
This VBS3/2 state can also be expressed as a matrix product 
states (MPS) 



IVBS3/2) = J2 Tr(S[*ilB['^l • • • B['"1)|zi 
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where } with m = 0, ±1, ±2 are given by the following 
4x4 matrices 
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m. CONTINUOUS QUANTUM PHASE TRANSITIONS 

In the model Hamiltonian Eq.(l), for J2 = 0, the ground 
state of the model is the VBSi state, [ ] while J3 = it cor- 
responds to the VBS3/2 state. [ ] When both J2 and J3 are 
non-zero, the model is no longer exactly solvable. Then, we 
expect that a quantum phase transition may be reached by ad- 
justing the value of J3/J2. To our knowledge, this is one 
of the few microscopic models, exhibiting continuous quan- 
tum phase transitions between two distinct topological or- 
dered phases. 



A. Ground state phase diagram 

In any one-dimensional quantum spin systems, the corre- 
sponding ground states can always be simulated by the wave 
functions in the MPS form, as the area law can be easily sat- 
isfied. To study the ground state properties of the general 
Hamiltonian Eq.(l), we thus propose a MPS with a finite local 
matrix dimension to approximate the ground state |V'g), 



Tr(. • • r"' Ar"'+iA • • • )| • • • m,m,+i • • • ; 



(12) 

where rrii = —2, —1, 0, 1, 2 is the spin quantum number of 
S'f , {r™'} is a set of X Z) dimensional matrices with bond 
dimension D corresponding to number of states kept in den- 
sity matrix renormalization group, and A is also a D x D 
dimensional non-negative diagonal matrix with its matrix ele- 
ments Xa, satisfying the normalization condition J^a — 1- 



The trace gives the superposition coefficients of the Hilbert 
space basis. The local matrices {F™} and A are set to be iden- 
tical on different sites due to the translation invariance.[20] 

The spirit of the infinite time evolving block decimation 
algorithm (iTEBDA) [lu] is to do the following evolution in 
imaginary time 



lim 



exp(-gT)l^o) 
\\cM-Ht)\^o) 



lim 



{exp{-He)f I*, 



(13) 

where j^E'o) is an arbitrary random initial state having nonzero 
overlap with \ijjg), r is imaginary time, and e is a small inter- 
val satisfying eN = t. It is a projection method with high ef- 
ficiency: as T increases, the weight of ground state in evolved 
state exp(— iJT)|5'o) grows exponentially. Next we split the 
Hamiltonian into two non-commutative parts. 



H 



+ + 



add 



z— even 



h{i,i + 1) — Hodd + Heven, 



(14) 



where all the local two-body bond operators h{i,i + 1) com- 
mute with one another in each part Hodd or Heven- Then 
Suzuki-Trotter decomposition is used to divide the respective 
evolution in sequential order, 

exp(-iZ£)|*o) = cxp{-eHodd) exp{-eHeven) + O(e^) 

(15) 

Owing to the fact that cxp{—eHodd) and cxp(— eiZeuen) com- 
pose of commutating evolution operators, all the bond evolu- 
tion in each part can be implemented simultaneously, which is 
compatible with the translation invariance of MPS. After each 
evolution step, the domain of local matrices increases by at 
least one site and the number of different local matrices in- 
creases by a factor at least 5. Then a singular value decompo- 
sition and bond dimension truncation are introduced to keep 
the MPS in the form given by Eq.(12) and bond dimension 
fixed in each evolution step, the detailed can be found in Ref. 
[16]. 

Now we show how to calculate the ground state energy den- 
sity. First we construct the transfer matrix 



G= J2 r"A®(r")*A 



(16) 



and compute its dominant eigenvalue 771, as well as the associ- 
ated right and left eigenvectors |ri) and \li). Then the matrix 
used to compute the energy expectation value is constructed 
as 

Ge= ^ (p,q|/i|s,t)rfAF«A® (r^)*A(r*)*A. (17) 

where (p, q\ and \s,t) are the wave functions defined by local 
Hilbert space of two adjacent sites and with the spin quantum 
numbers p, q,s,t = —2, —1, 0, 1, 2. The ground state energy 
per site Eg can be calculated by 



Eg = lim 

" JV-i-oo 



tr (G^) 



(18) 
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FIG. 1 : (Color online) Ground state energy density varies as a func- 
tion of Jg for fixed J2. The local matrix dimension D is set to be 
300 in the calculation. 



FIG. 3: (Color online) The second order derivative of the ground 
state energy density varies as a function of J3 for fixed J2 . 
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FIG. 2: (Color online) First derivative of the ground state energy 
density varies as a function of J3 for fixed J2 . 



FIG. 4: Ground state phase diagram of the model Hamiltonian 
Eq.(l). 



For a fixed J2, we calculate Eg as a function of J3. The 
quantity Eg and its first order derivative with respect to J3 
are finite and continuous, which are displayed in FIG.l and 
FIG. 2. However, the second derivative of Eg with respect 
to J3 exhibits divergence as J3 is tuned to a certain critical 
value. This is similar to the specific heat divergence in classi- 
cal phase transitions. Such a behavior is shown in FIG. 3 for 
several typical values of J2. We thus conclude that the system 
undergoes a second-order phase transition at zero tempera- 
ture. By determining the positions of the critical points, we 
thus derive the zero temperature phase diagram in FIG. 4. 

Moreover, the calculations of the spin-spin correlation 
length and entanglement entropy also show a singular behav- 
ior and provide further evidence of the second order phase 
transition of the model. In FIG. 5, the numerical results of both 
spin-spin correlation length and entanglement entropy are de- 
picted as a function of J3 for fixed J 2- At the critical point, the 
spin-spin correlation length is divergent, and the entanglemnt 
entropy shows a cusp. Away from the critical point, an extrap- 
olation of correlation length and entanglement entropy can be 



calculated as D goes to infinity, indicating that both of them 
saturate to finite values. Due to the finite spin-spin correlation 
length and the form of the model Hamiltonian in terms of pro- 
jection operators, it is straightforward to prove the existence 
of excitation gap[2 1 ] and to identify two phases separated by 
the critical line in the J2 — J3 phase diagram as VBSi and 
VBS3/2 states, respectively. The computational methods of 
entanglement entropy and correlation length are explained in 
detail in the following sections. 



B. Entanglement spectrum across the transition 

Li and Haldane[3] have recently proposed that entangle- 
ment spectrum (ES), i.e., the minus logarithms of the eigen- 
values of a reduced density matrix, can be used to characterize 
topological order If there is an entanglement gap separating 
the low-lying ES and the upper parts, then one can find a one- 
to-one correspondence between the low-lying ES and the low 
energy spectrum of individual edge excitations. In particu- 
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FIG. 5: (Color online) (a) The spin-spin correlation length varies as 
a function of J3 for two typical values of J2 . (b) The entanglement 
entropy varies as a function of J3 for two typical values of J2. The 
local matrix dimension D is fixed at 300 in this calculation. 

lar, the lowest level of the entanglement spectrum for a topo- 
logical ordered state should be degenerate. When the state 
changes from VBS3/2 to VBSi by tuning the ratio of J3/J2, 
how does the topological order changes in this process, espe- 
cially when crossing the phase transition point? We will try to 
use entanglement spectrum as a probe to partially answer this 
question. 

If the MPS in Eq.(12) is in "canonical form", then upon 
dividing the system into left and right parts the ground state 
wave function should become 

|vI/,)=^A„|$^)|$f), (19) 

a 

where {|$^),« = 1,2,---,D} and = 
1, 2, • • • , D} are orthogonal basis states of the left and right 
semi-infinite chain, 

($^|$^) = (<f>fici>^')==(5„;3 (20) 

It can be shown that the canonical condition Eq.(20) imposes 
the following constraint[ ] on the F'" and A in Eq.(12), 

^r™A2(r")^ ^ ^(r™)tA2r" = Id^^d^- (21) 

In general, a MPS has a gauge freedom, because the local 
matrixes can be the same up to a similarity transformation. 



FIG. 6: (Color online) The thirty-sixth lowest values of the entan- 
glement spectra for the VBS 3/ 2 and VBSi phases on two sides of 
the critical points. (a)(J2,J3) = (0.1, 0.029), (b) (J2, J3) = 
(1.0,0.25), (c) (J2, J3) = (3.0,0.53). A light gray vertical line 
is put on the critical point to separate the two phases. Each entan- 
glement spectrum is represented by a small cross, and the degenerate 
spectra are spatially staggered a little bit in horizontal direction to 
distinguish and count them. In this calculations, the local matrix di- 
mension D is set to be 300. 



In canonical form, it is extremely easy and straightforward to 
write down the entanglement spectra Pa = — log(A^) . In 
our calculation, we perform the canonical transformation[23] 
explicitly at the end of iTEBDA and obtain the MPS in its 
canonical form. 

The 36 lowest values of ES are plotted in FIG. 6 for J2 = 
0.1, 1.0 and 3.0. It can be clearly seen that, for the VBS3/2 
state, the lowest entanglement eigenvalue is four-fold degen- 
erate in one to one correspondence with the four-fold de- 
generate edge states; while for the VBSi state, the lowest 
eigenvalue of the ES is three-fold degenerate. Above this de- 
generate eigenvalues there exists a large gap, so the degen- 
erate levels are protected topologically. So such a calcula- 
tion of the ES can be used to confirm that both VBSi and 
VBS3/2 are really topological ordered states, as well as that 
the ground state degeneracy of a long enough open chain are 
really 16 and 9 respectively. When approaching the critical 
points, the degeneracies are gradually lifted but the gap still 
survives. However, in the vicinity of quantum critical points 
(J2, J3) = (0.1,0.029), (1.0,0.25), (3.0, 0.54), the degenera- 
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cies in the ES no longer exist and the large gaps between the 
degenerate lowest level and the higher levels are no longer 
present. Due to the expected finite-size effect near the critical 
region, so far we can not simply conclude that the topolog- 
ical order is destroyed completely at the critical points. Re- 
cently a partition with a very non-local real space cut has been 
proposed[24], and some new light has shed on using ES to de- 
tect non-local orders in gapless spin chains. However, it still 
needs further investigation to clarify this question and will be 
done in future works. 



C. Central charge on the critical line 

At the critical points, the system should be described by 
conformal invariant quantum field theories. For such theories 
the central charge encodes information about the universality 
class. According to conformal field theory, the von Neumann 



entanglement entropy should diverge logarithmically with the 
correlation length[2?], 

Se^^HO + So, (22) 

where Se is the entanglement entropy between two semi- 
infinite parts of a whole chain, c is the central charge, ^ is 
spin-spin correlation length in units of lattice spacing, and So 
is a non-universal constant. For a MPS in canonical form, 
entanglement entropy can be calculated easily[ 14, 15]: 

a 

where are the coefficients in Eq.(19). 

The spin correlation length can be deduced by two points 
spin-spin correlation function as 
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where transfer matrix G, eigenvalues 771, and eigenvectors {li \ 
and |ri) have the same definition as in Eq.(16). 772 is the sec- 
ond largest magnitude eigenvalue of G and (^2! and \r2) are 
corresponding left and right eigenvectors. Similar to Gb in 
Eq.(17), Gz is defined by 



G, = ^(p|^^|<7)rfA® (r«)*A. 



(25) 



Generally, for a state without spin long-range order, 
{li\Gz\ri) ~ 0, the two points correlation function can be 
written as 



lim lim {S^S' 



(26) 



where the correlation length ^ = 1/ log(| 771 7772 1). 

Using Eq.(23) we have calculated Se and ^ in the vicinity 
of several different critical points on the phase boundary of 
FIG. 4. The associated scaling relation between Se and f are 
shown in FIG. 7. Although there may be some deviations at 
small ^, 5*,. tends to lie on the c = 2 line for large ^. The fact 
that all the central charges are approximately equal to 2 im- 
plies a single fixed point governing the critical behavior of the 
entire phase transition line. So far the conformal field theory 
with c > 1 can not be classified systematically, and there- 
fore to determine the corresponding conformal field theory of 
a fixed line with c = 2 might be worth attempting. 

Here we present some conjectures deduced from the con- 
formal field theory kinematics. According to the central 



charge value and the constituents of the VBSi and VBS3/2 
states, the level-four SU(2) Wess-Zumino-Witten model is the 
most possible effective field theory and the conformal weight 
of the primary field given by[26] A*-^-* = j{j + l)/6 with 
j = 0, 1/2, 1,3/2,2. Actually, the level-four SU{2) Wess- 
Zumino-Witten model can be regarded as the effective critical 
field theory of the following spin S ~ 2 antiferromagnetic 
Takhtaj an-B abuj ian model [ ; ] : 



-^^(SjSi+i)^ - — (SiSi+i)'' , (27) 



which can be written in terms of the projection operators as 

H = J^[JiPi(i,i + l) + J2P2(j,« + l) 
i 

+ J3P3{l,i + l) + P4{l,l + l) (28) 



22 
25' 



with >/i = M, ^2 = § and J3 

Compared to the model Hamiltonian Eq.(l), there appears 
an additional interaction term Pi{i,i + 1) with the largest in- 
teraction strength Ji = i|. By calculating the entanglement 
spectrum, we find that this critical point just lies on the bound- 
ary between the AKLT and a dimerization phases. For the 
dimerization phase, the entanglement spectrum show an even- 
odd difference from the topological ordered phase, i.e., if the 
bipartition is done at a bond connecting left even and right odd 
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FIG. 7: (Color online) The scaling behavior of the entanglement 
entropy 5e as the correlation length ^ in the vicinity of different 
critical points, (a) (0.1,0.029), (b) (1.0,0.25), (c) (2.0,0.42), (d) 
(3.0,0.54). The black solid line is as a reference with a central 
charge c = 2 and the red squares are numerical results. In our calcu- 
lation, the local matrix dimension increases from 12 to 600. 

sites, the lowest spectrum is 5 fold degenerate, otherwise the 
lowest spectrum is non-degenerate and has a big gap with the 
upper part. This even-odd difference indicates the existence 
of dimerization phase. The relation between this critical point 
and the 50(5)-AKLT critical line can be understood as fol- 
lows. When we fix Ji = il and J2 = ||, J3 varies from 
to ||, the system evolves from the 50(5) symmetric phase to 
the dimerization phase, and then the AKLT phase for J3 > ||. 
Moreover, for the fixed value of J2 the dimerization region 
shrinks and finally disappears when Ji is decreased, which is 
compatible with our ground state phase diagram. Therefore, 
we expect that there exists a crossover flow from the fixed 
line of the transition between the 50 (5)- AKLT phases to the 
S — 2 antiferromagne tic Takhtajan-Babujian model. 



IV. CONCLUSION 



In summary, we propose a one-dimensional spin-2 Hamil- 
tonian, which exhibits two topologically distinct VBS states 
in different solvable limits. By using the infinite time evolv- 
ing block decimation algorithms, we have studied the quan- 
tum phase transition between them and determined the central 
charge to be c = 2. Of course, continuous phase transition be- 
tween topological phases characterized by different number of 
edge states is known. For example, by tuning the coefficient of 
the topological term in the 50(3) / S0{2) non-linear a model, 
it is possible to induce phase transition between VBS states 
associated with different spin values. The transition studied in 
this paper is very different. It takes place between two topo- 
logically distinct VBS states associated with the same spin 
value. We are not aware of any previous study of this type of 
phase transition. 
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Note Added After we submitted the original version of this 
manuscript for publication, a paper[28] concerning with the 
similar issue by different numerical density matrix renormal- 
ization group method on a finite length of chain appeared on 
the archive, where the authors claimed the existence of dimer- 
ization phase separated the VBSi and VBS3/2 phases in the 
ground state phase diagram. However, if the ground state 
energy density and its second-order derivative do not show 
any singularity as the coupling parameters approach to the 
boundary of the "dimerization phase" on both sides, no quan- 
tum phase transition can occur, and the claim of dimerization 
phase existence is not reliable. 
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